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Abstract 

Based on the multiple scattering series incorporated in the Faddeev scheme 
the high energy limit of the total n+d cross section is evaluated in a non- 
relativistic model system where spins are neglected. In contrast to the naive 
expectation that the total n+d scattering cross section is the sum of two 
NN cross sections we find two additional effects resulting from rescattering 
processes. These additional terms have different signs (shadowing and anti- 
shadowing) and a different behavior as function of the energy. Our derivation 
of these results which are already known from Glauber theory is based on 
the analytical evaluation of elastic transition amplitudes in the high energy 
limit. It does not depend on the diffraction type assumptions connected with 
Glauber theory. In this model of spinless Yukawa type forces (with no ab- 
sorption) the total n+d cross section does not approach twice the NN total 
cross section in the high energy limit but rather approaches the total NN cross 
section multiplied by a number larger than two. Therefore, the enhancement 
effect resulting from rescattering is larger than the shadowing effect, which 
decreases faster with energy. 
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I. INTRODUCTION 



Recently it became possible to calculate the total cross section for neutron-deuteron 
(n+d) scattering with high precision in the energy regime up to 300 MeV projectile en- 
ergy by solving the nonrelativistic 3N Faddeev equations based on modern NN forces [|l|. 
Compared to the most naive picture, which in the high energy limit would equate the to- 
tal n+d cross section with the sum of the total cross sections for neutron-proton (np) and 
neutron-neutron (nn) scattering, the rigorously calculated result up to 300 MeV is smaller. 
Obviously one can expect some shadowing effect in the reaction, which would explain this 
result. On the other hand, rescattering of the nucleons upon each other might a priori 
also enhance the total n+d cross section over the sum of the individual two nucleon cross 
sections, especially if the forces are attractive. In the context of Glauber theory [^]-|^ both 
features, enhancement and weakening, are present. In a model of spinless particles, which is 
based on the Faddeev formulation, we want to study the high energy limit of the total n+d 
cross section and evaluate the leading terms analytically. This approach differs from the one 
used in the Glauber formulation. Based on the multiple scattering series incorporated in the 
Faddeev framework, we calculate the first and second order terms of the series when taken 
in the high energy limit and show the contributions of the two terms to the total n+d cross 
section. Starting from a multiple scattering series implies that this is an ordering according 
to powers in the NN t-matrix. Though we restrict ourselves to a nonrelativistic framework, 
we nevertheless consider the results as being instructive. The analytical steps leading to 
the high energy limit, which only involve ordinary analysis, are carried out in well defined 
integrals. There are no a priori assumptions about the scattering process involved, like e.g., 
diffraction type approximations. We also take the identity of the particles explicitly into 
account. In this paper the complications due to the inclusion of the spin and isospin degrees 
of freedom are avoided so that the basic mechanism can be seen more clearly. Taking spin 
and isospin degrees of freedom into account will modify the results due to the interferences 
of spin and isospin dependent terms, as will be shown in a forthcoming article. Specifically 
because of those spin and isospin interference effects, which are quite involved, we want to 
present this more transparent case with three bosons separately. 

In Section II we describe the Faddeev framework, its multiple scattering expansion, and 
the leading order terms in the NN t-matrix for obtaining the total n+d cross section. The 
high energy limit of the corresponding expressions is carried out analytically in Section III. 
To illustrate the behavior of the leading order terms in the high energy limit numerical 
examples are given in Section IV for a superposition of Yukawa interactions. We conclude 
with Section V. 



II. LEADING MULTIPLE SCATTERING TERMS FOR THE TOTAL N+D 

CROSS SECTION 

We consider three identical spinless bosons which interact via two-body forces. In our 
usual manner p to exploit the Faddeev scheme the operator for elastic scattering of a 
nucleon from a bound nucleon pair is given by 

U = PGo^ + PT, (2.1) 
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where the three-body operator T obeys the Faddeev type integral equation 



T|$) = tP|$) +tPG'oT|$). (2.2) 

The channel state, which is composed of a two nucleon bound state (deuteron) and a mo- 
mentum eigenstate of the projectile nucleon, is denoted by Furthermore, t is the two 
nucleon transition operator. Go the free three nucleon propagator, and P the permutation 
operator, which is the sum of a cyclic and an anticyclic permutation of the three particles. 
The elastic forward scattering amplitude is given by the matrix element ($|?7|$). When 
iterating Eq. (|2.2| ) and inserting the result into Eq. we obtain the multiple scattering 

series for the elastic forward scattering amplitude 

($|[/|<|.) = ($|PGo + ($|PtP|$) + {^PtPGQtP\^ 

+ {^\PtPGotPGotP\^) + ■ ■ ■ , (2.3) 

which is an expansion in orders of the NN t-operator. Using the optical theorem one obtains 
for the total cross section for nucleon-deuteron (n-|-d) scattering 

4t77 
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Here Qq is the asymptotic momentum of the projectile nucleon relative to the bound two-body 
subsystem. From Eq. (2.3) follows 

2i Jm($|f/|$) = ($|f/|$) - ($|f/|$)* 

= (<l>|P(t - t^)P + PtPGotP - Pt^GlPt^P + PtPGotPGotP 

-Pt^G*oPt^G*oPt^P\^) + ■■■ . (2.5) 

Since the first term in Eq. (2.3) is real, it does not contribute to the total cross section. 

For the analytical extraction of the high energy limit it is useful to rewrite Eq. (2.5) in 
the following form 

(<l>|f/|$) - ($|f/|<l>)* 

= ($|P(t - t^)P|$) + ($|P(t - t^)PGotP\^) 

t - t^)PGotP\^y + ($|Pt^P(Go - G*)tP\^) 
t-t^)PGotPGotP\^) 
^P{Go-G*)tPGotP\<^) 
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($ 



^PG*t^PG*{t-t^)P\<!)) + ■■■ . (2.6) 

The next step is to explicitly evaluate the permutations given by P = P12P23 + -^13-^23- This 
specific choice of the permutation operator P corresponds to the choice of particles 2 and 3 
forming the two-body bound state and 1 being the projectile, i.e., t = ^23 and 

1$) = |$)i = |vP,)23|qo)i. (2.7) 
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The subscripts denote which particles occupy the states. This specific choice leads to 

|$)2 = ^12^231$)! = |V'd)l3|qo)2 (2.8) 

and 

|*)3 = Pl3P23\^)l = |^d)l2|qo)3- (2.9) 

Taking advantage of the symmetry property of the subsystem, P23|'^')i = l*^')! and P23^23 = 
^23-^23, one obtains after some algebra 

($|P(t-tt)P|$)=2 2($|(t-t^)(l + P23)|*)2 (2.10) 

and 

($|PtPGotP|$) = 2 2mtGot2 (|$)i + |$)3) + 2 2mtGoh (|$)i + 1^)2) . (2.11) 

A similar evaluation could be done for the terms third order in t given in Eq. ( ^.5] ). For 
our present considerations of the n+d total cross section in the high energy limit we restrict 
ourselves to the terms given in Eqs. ( 2.101 ) and (|2.11 ). 



III. THE HIGH ENERGY LIMIT 



First we need to consider the exact momentum space representations of Eqs. ( 2.10|) and 
( [^.11| ). Equivalent decompositions of the unity operator are given by 

1 = J £p d^q \p(l)i i{pq\ , i = 1,2, or 3. (3.1) 

Here |pq)i = |p)i|q)i and p and q are the standard three types of Jacobi momenta for three 
particles 0. The index i denotes the singled out nucleon. Inserting the unity operator into 
the first term in Eq. ( ^.lOj ) leads to 



2($|i-t^|$)2 = / d^pd'q J d^p'd'q'J dYd'q" J dY'd'q'" 2(<f|pq)2 2(pq|p'q)i 

X ,(p'q'|t - tt|p"q"), i(p"q"|p'V')2 2(p'"q'"|$)2. (3.2) 

Further, one has 

,(pq|<l>), = (5=^(q-qo)(p|¥^d). (3.3) 
The standard relations among the different sets of Jacobi momenta give 

i(pq|p'q')2 = S^p - iq - q') 6'{p' + q + iq'). (3.4) 



In addition one has 



L(pq|t(i?)|p'q')i = {p\tiE - £g^)|p') - q'). (3.5) 
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Employing all the above given relations, it is straightforward to arrive at the following 
expression for Eq. (3.2) 

2($|t-t^|$)2 = J d^q (v^dl -q - |qo) (|q + qo|t-t^||q + qo) (-q- Iqclf^d) 

= J d'k {^,\ - k) (-k|v,,) (fqo + ik|(t - tt)(5)||qo + ik). (3.6) 



The total energy is E = ed + ^q^, where is the deuteron binding energy. When expressed 
in terms of k = q + ^qo one obtains for the energy argument e of the t-matrices 

^ . 9 /I 1 

2^0) 



^ = -^- I^^ - ed + j^qo - 

2 



Q+^(fqo + |k) -^k 



Am ^ 
„2 



(3.7) 



If the projectile momentum go is sufficiently large in comparison to the dominant deuteron 
momenta contributing to the integral in Eq. (3.6) and thus e ~ ^(|qo + |k)^, one encounters 
on shell NN forward scattering amplitudes under the integral. The permutation operator 
P23 in Eq. ( glOD leads to the necessary symmetrization of the t-matrix elements 



where 



|q)s = |q) + I - q), 



(3.8) 



(3.9) 



and q = |qo + |k. Using the Lippmann-Schwinger equation for t, with the two-body force 
V as driving term, and q = |g|q one has 

(q \{t-t^){e=^q^M, = -2vrz^g J dq' (q|l^|gqV+^ (gq1^|q)« 

= -27ii — q J dq {q\t\qq) s{q\t\qc^)* 
= -ni^q J dq' |(q|^|gq%|^ 

This expression is directly related to the two-body total cross section 



(3.10) 



NN 
tot 



(q|ikq')s^ 



(3.11) 



We adopted here the usual convention for the total cross section for two identical particles 
IHl and obtain 



- mis = ^q'M 



-iq 



1 



NN 



m (27r^ 



3 ^tot 



(3.12) 



As a consequence Eq. ( |3.6| ) supplemented by the symmetrization as given in Eq. ( |2.1CI| ) and 
approximated in the energy argument of the two-body t-matrix can be rewritten as 
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2($|(t-t^)(l + P23)|$)2 (3.13) 

If the projectile momentum go is sufficiently large compared to the typical deuteron momenta 
the function ||qo + |k|cr^/^(^(|qo + |k)^) is expected to vary slowly over the range of the 
k- values contributing to the integral in Eq. (3.13). Using the normalization of the two-body 
bound state, expanding the function at k = and knowing that the contribution of the first 
derivative vanishes, one obtains in the limit for large qo 



im - tt)(i + p,sm, -^j^sl ^o^iTf (^(i^o)^) / (^,1 - k) (-k|^,) 

' T^oaiTf (^(fgo)^). (3.14) 



m (27r) 
2i 1 
~m(2^34 

The arguments leading to Eq. (3.14) are similar to those used in arriving at the method of 
optimum factorization successfully applied in intermediate energy pion-nucleus and nucleon- 
nucleus scattering P,p!0|. The n+d total cross section, Eq. ( p.4| ), thus gives in the first order 
term in t of Eq. (2.6) in the high energy limit 

(^tot 1 1st order ^ 2 cr^/^. (3. 15) 

This result corresponds to the naive expectation that at high energies the projectile 
nucleon sees the individual nucleons inside the deuteron as if they were independent particles. 
It should be pointed out that due to the optical theorem a^^l^ is O(t^), though the expression 
of Eq. (|3.15|) has been derived from the terms linear in t in Eq. ( |2.6|) . 

Let us now consider the contributions to the total n+d cross section second order in t as 
given in Eq. (p75|). A straightforward but somewhat tedious algebra using Eq. (2.11) leads 
to the exact form 

{^iPrPGotPl^) 

{<^d\q}{q.'\Vd) 



i 9'^ + q-q') + 2^(q + q')-qo + «e 

xdqo + iq|r(ei)||qo - h " q') (-fqo + q + iq1^(e3)||qo + ^q% 

{fdh){(l'\<fd) 



+ 2 Jd^qJ d^q' 



IqI - ^(g^ + g'^ + q ■ q') + al^ (q + q') ■ qo + 

X (f qo + |q|r(ei)| - fqo + iq + q') (f qo - q - |q'|t(e2)| - f qo - |q% (3.16) 

The quantity r represents either t — or t"^ , and thus acts in the subsystem 1=(23). Again, 
the subscript s indicates the symmetrized state as given in Eq. (|3.9|) . The energy arguments 
in the t-matrices are 

^^ = ^d + i (fqo)' + il^q ■ qo - (3.i7) 

and 

e3 = e, = e, + ^ (fqo)' + £q' ■ qo - £q" (3.18) 
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It should be pointed out that the occurrence of the symmetrized state incorporates the 
forward and backward scattering amphtudes. 

For the hmit ^ oo one can again neglect the variations of the t-matrices under the 
integrals and obtains in this limit 



{^iPrPGotPm ^ 2 (|qo|r(i^(|go)^)l!qo). (|qo|t(^(|go)^)l!qo). x J, (3.19) 



where 



f d\' , , {^M){^W.) ^ _ (320) 

-|ed|-^(g2 + g'2 + q.q.) + _^(q + q/).q^ + ,^ ^ > 

The interesting insight will now arise from the quantity J, which we still need to consider in 
the high energy limit. These considerations are most transparent if one works in configura- 
tion space. The momentum space form will be given in Appendix A. Let 

>,<\Wd) = jtA^ [ d\ e'^-' (rl^rf) (3.21) 



with 



(27r)3/2 



(r|^d) = -^^d(r). (3.22) 



Then the quantity / from Eq. (|3.20|) becomes 



(2vr) 

-kdl - ^(g' + + q ■ q') + 2l^(q + q') ■ qo + 

Since the denominator depends on qo, the quantity whose infinite limit we want to consider, 
it is natural to split the vector q (q') in components parallel and perpendicular to qo, 
and consider qo pointing into the z-direction. For the vector r (r') we define a similar 
decomposition 

q=(q±,g.) (3.24) 
r = {rj_,z). 

and obtain 

/ = / d'q^ J d'q'^ J d'r y,,(r) J d'r' ^,{r') e^^-'^- e-'^'-""'- (3.25) 



(27r)3 47r 
X J dqz J dq^' e*^""' e" 



-m\ed\ - gi - q'l - q± ■ q'± - - q'l - QzQz + lqo{qz + q'-,) + ie' 
Further we substitute 
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qz + q, = s 



(3.26) 



so that the denominator takes the form 



with 



= m\ed\ + q± + q''± + q± ■ q'± + s'^ > 



/2 



(3.27) 



(3.28) 



For large go one obtains 



D ^ -|(s - 2qo - is){s + is) 

<j qo 

-|(s - 2qo - ie){s + ie) 



(3.29) 



and thus encounters two poles in the complex s-plane. The integrals over s and s' can be 
carried out analytically, leading to 



oo oo 



ds / ds' e^(^+^')-'ef (--^')- 



2i 



-4/3 



-oo — oo 



's — 2qo — ie){s + ie) 



- — (27r)^5(z + z') e(z)e^'''°' + e(-z) 
3go '- 



(3.30) 



Inserting this into Eq. (3.25) and carrying out the integrations over and q^ yields 

oo oo 

I^m (27r)2 ^ J d\± J d^r\ j dz j dz' (3.31) 

— oo — oo 

X <pa + 5{v^) 5{v\) 5{z + z') [e{z)e'''i°^ + e{-z) 



1 

47r' 



The (5-functions under the integral show that contributions to / only arise if the two nucleons 
in the deuteron sit behind each other with respect to the projectile momentum qo = qoz. 
This coincides with our naive understanding of shadowing. The term containing e^'^°^ falls 
off fastest in Eq. (3.31) and will be neglected. Thus, one arrives at 



1 °° 1 

^ 

-(27r)3((^,|-|(^,) — . 



3go 



(3.32) 



The corresponding algebraic steps performed in momentum space are shown in Appendix A 
and yield the result 



(3.33) 
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where p> = max{p,p'). 

We now return to Eq. (3.19), set r = t — t^, use the relation given in Eq. ( p.l2|) to the 
total NN cross section and obtain 

($|P(t-tt)PGotP|$) - -2airf (|qo|t(^(fgo)^)||qo). (v^.l^lv^d)^- (3-34) 

Finally subtracting the conjugate complex according to Eq. (2.6) and using the optical 
theorem in the two-body subsystem, 

(fqoltlfqo). - (fqoltlfqo): = -2^^ {^)' <^ ^ (3-35) 



($|P(t - t^)PGotP\<^) - ($|P(t - t^)PGotP\<^y 



' i-fon' (3-36) 



leads to 



Because of (cr^/^)^ this expression is of order 0(i:^). 

It remains now to discuss the last term of second order in t in Eq. (2.6) in the high 
energy limit. Here we encounter Go — = —2m5{E — Hq) and obtain 

($|Pttp(Go-GS)tP|$) 



(!qo|t(^(!go)')||qo).|7 d'q I d\' (^,|q) (q>,) 



X 5 - ^(g' + + q ■ q) + olrqo ■ (q + q')) • (3-37) 



The integral term can simply be related to the expression I from Eq. (3.20). Similarly to 
Appendix A we work in momentum space and obtain 

I'^j d'q I d'q' (^,|q) (q^,) 5 {-\e,\ - l.(q^ + q'^ + ^ ■ q') + ^qo ■ (q + q')) 

= / d^q± J d\'^ J dq,J dq,, {^Jql + ql) Vd {^^q'l + q'l^ ^ 

X 6 {-\e,\ - ^{ql + q'l + q^ ■ q'^ + q^ + q']) + ^^loi^z + q'S) ■ (3.38) 
Then using Eq. (3.29) for the argument of the 5-function we arrive at the dominant term 

I'^mJ rfV / d\'^ J ds'J ds^^Jis)^, ^ql + (.' + i.)^) ^q'l + (s' - i.)^) i- 

= ^ / d'q^ J d\^ J ds' (vgiTT^) V, (vViT^) ^ 

1 

/. (3.39) 



The last step follows from Eq. (A.l), if one evaluates this integral in a different way. Thus 
we end up with 
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($|Pttp(Go-GS)tP|<|.)^-2z — (2vr)3 (3q^|t(i.(3^„)2) 1 3q^)^ (3.40) 



2 , , 1 , , 1 



Adding the above expression to Eq. ( |3.36|) and using Eq (2.4) gives the asymptotic contri- 



bution of the terms second order in t of Eq. (2.6) to the total n+d cross section as 

1 2nd order = -2 (cr^^f) (v^d | — | '/'d) 2" (3-41) 



'a \ 2 

4m \ 



3go/ 



+ (2vr)« ( ^ ) (|qo|t(^(|go)')l!qo)s ' (^dl^l^d) 



1 , > 1 



7-2 4^ 



If we use the optical theorem in the 2-body system for the imaginary part occurring in 
l(!qo|t(^(|go)^)|fqo)«P, we obtain 



(^tot\2nd order ^ " (o"tof ) (Vdl^l'Pd)^ 

2 



+ (2^)' (^) (^e(|qo|t(i,(|go)^)l!qo).)' (^.I^lv^^)^- (3-42) 

In summary, we obtain in the high energy limit for the total n+d cross section based on 
the terms of first and second order in t in the expression for /m(<l>|t/|$) 

'^tot ~ O toil 1st order "T Ojof 1 2nd order 



NN 

2 



+ (2-)^ (^) (i?e(fqo|t(i(fgo)^)||qo).)' (v'^l^^)^ 

= 2aZ^ + 0it')-0{t'). (3.43) 

As mentioned earlier, the first term, obtained from the first order term in the two-nucleon 
t-matrix within the multiple scattering expansion gives twice the NN total cross section, 
thus considering the two nucleons in the deuteron as being free particles. The terms of 
second order in t within the multiple scattering expansion of Eq. ( |2.6| ) give rise to two 
correction terms. The positive term in Eq. (3.43), being of O(t^), enhances the contribution 
of the first term. It is tempting to consider this term as an anti-shadowing effect. The 
third, negative term in Eq. (3.43), which is due to (cr^/^)^ of Oit"^), reduces the value of 
the two total NN cross sections and its action is naturally called a shadowing effect. In the 
Glauber approximation one arrives formally at the same result. However, it should 
be pointed out that Eq. (3.43) contains the symmetrized two-body t-matrix elements and 
thus implicitly forward and backward amplitudes. We think that our derivation based on 
the multiple scattering expansion of the Faddeev formulation together with the momentum 
space treatment is a viable alternative. It uses standard and transparent evaluations of 
integrals, which could also be calculated numerically without going to the high energy limit. 
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IV. APPLICATION 



In this Section we want to give a numerical illustration of the results derived in the 



previous section. We use two NN model forces of Malfiiet-Tjon type |Tl|, one being purely 
attractive and the other having attractive and repulsive parts. Both are of Yukawa type and 
given as 

The parameters are given in Table I. As we saw in the last section, there are two momenta 
which control the asymptotic expansion. The first is the initial projectile momentum qo, 
which appears as |qo in the three-body context. The second is the 'typical' deuteron 
momentum. For the two model potentials we display in Fig. 1 the corresponding deuteron 
wave functions. We see that in both cases the wave functions drop by about a factor of 
10 within 0.6 fm~^. Thus we expect that for |go sufficiently larger than this value the 
asymptotic expressions derived in the previous section should be vahd. In order to get 
quantitative insight into the onset of the validity of the asymptotic expressions one should 
evaluate the terms in the multiple scattering series exactly. Strictly spoken this amounts to 
solving the Faddeev equation, Eq. (2.2). As a first step we exactly evaluate here the term 
in Eq. (2.10), which is in first order in t. Using Eq. (3.6) one finds the following exact form 

(<l>|P(t - tt)P|$) = ^ [d'p ^^(2|p - fqol) Im{p\t{e)\p)s (4.2) 

where 

e = e, + ip'-^{p-lqor (4-3) 

Introducing as variables explicitly the magnitudes of the momentum vectors p and qo 
and the angle between them, x = p ■ qo, we obtain for the n+d total cross section in the 
first order in t 

4 1 oo 

^rilistorder = -(277)^^8 J dx J dp (^2^p^ + (|go)2 " Ipqox^ Imt,{p,p,e). (4.4) 

-1 

The integral over the imaginary part of the forward scattering amplitude receives contri- 
butions from the region e > and the deuteron pole. The condition e > limits the 
integration regions in x as well as p. The x-integration has to be carried out only between 

Xmin = Jl + and ^rnax = 1, and p-values given by 



Pmax,min = QqX ± \/ " f^O ^ 3 I I • (4-5) 

The pole contribution is obtained from an integration over x with similar boundaries at fixed 
p- values. 

The off-shell NN scattering amplitude is determined by firstly solving the two-body 
Lippmann-Schwinger equation exactly in three dimensions, i.e. without partial wave expan- 
sion ||12| . The quantities entering the integral in Eq. ( [4.4|) are off-shell t-matrix elements for 
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equal magnitudes of the momenta. Thus the integration requires a two-dimensional inter- 
polation of the imaginary part of the forward scattering amplitude in the variables p and e. 
We use standard b-splines for the numerical calculation . 

The first order term of the total n-|-d cross sections for the two potential models MT-III 
and MT-IV is shown in Figs. 2 and 3 as function of |go and compared to its asymptotic limit 
2c"tot^ given in Eq. (3.15). A closer inspection reveals that this limit is reached at nucleon 
laboratory energies of about 300 (600) MeV within 1% for MT-III (MT-IV). This energy is 
fairly high so that the simple potential picture will no longer be valid and relativistic features 
will be important, including meson production. These absorption processes will of course 
change our results as is known from studies within the Glauber formalism . Nevertheless 
we think it should be of theoretical interest to see the results within a pure potential picture. 
We leave the exact solution of the Faddeev equation, Eq. (2.2), to a future study and estimate 
now the resulting rescattering terms via the asymptotic expressions contained in Eq. (3.43). 
They are expected to be a good approximation around 300 MeV and above for the potential 
containing the realistic feature of repulsion. 

The deuteron matrix elements for the two potentials entering in the asymptotic expres- 
sion of Eqs. (3.32) and ( p.43| ) are given in Table I. 



We start our investigation with the MT-III potential, which has a short range repulsion 
and a intermediate range attraction, and thus contains some realistic features of the NN force. 
In Fig. 4 the three different terms of the right hand side of Eq. (|3.43|) are shown separately 
as function of the nucleon laboratory energy together with their partial sums building up 
the n+d total cross section in the high energy limit. We see that the rescattering term of 
O(t^) is quite large and can not be neglected in relation to the leading term 2a^^f in the 
whole energy range shown. For a more quantitative inspection the three different terms are 
listed separately in Table II for the higher energies. The first column gives twice the NN 
total cross section, the second column shows the positive term, which is also of 0{t^) like 
the total NN cross section and the third column gives the negative term, which is of O(t^), 
and which causes shadowing. We see from the table that the positive term of O(t^), which 
results from the second order in t of the multiple scattering expansion of the elastic forward 
scattering amplitude is relatively large. Already at 100 MeV the shadowing term 

is much smaller than the other two terms and drops of course quickly due to its energy 
dependence O(^). From these numbers we see that the rescattering corrections are quite 
large in this pure potential picture. This is even more pronounced for the purely attractive 
potential MT-IV as shown in Fig. 5. The corresponding terms are displayed in Table III. We 
see that the positive term of O(t^), which results from the second order in t of the multiple 
scattering expansion of the elastic forward scattering amplitude ($|?7|$) is quite large and 
more important than the sum of the two NN total cross sections. 

In the purely attractive potential model MT-IV we are able to rigorously provide the 
extreme high energy limit of the n-|-d total cross section. In this limit the Born term for 
the NN t-matrix gives the only contribution and we can replace the t-matrix element in 
Eq. (3.41) with the potential matrix element, which is real, 

(!qo|t(^(!go)^)||qo)s - (|qo|V^||qo)s. (4.6) 

This expression can be related to the total cross section, and one finds after a simple inte- 
gration 
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/ #KigoP|V^liqo)s|'^47rA^^(3q^|v|3q^)2. (4.7) 



Using the definition of tlie NN total cross section we obtain 



rqoM^qoL-l-j (^^--- (4. 



As a consequence one can write Eq. ( |3.43|) as 

[-^oH^ iv^l^M)^- (4-9) 



For tlie potential MT-IV the term additive to 1 is 3.166 (using the quantities given in Table I) 
and therefore the total NN cross section is multiplied by 8.33 if the exact limit is reached. 
This number is substantially larger than the naive expectation, which would be twice the 
NN cross section. Numerically this asymptotic limit is reached around 4000 MeV nucleon 
laboratory energy within 2% as we see from the last column in Table III. There we display 
the ratio of the high energy limit obtained from Eq. ( |3.43| ) and the exact limit from Eq. ( |4.9| ) 



R 



2ailf +(27r)6 


\ 3go ) 


' (^e(|qo 


^ m 


(!?o)^)||qo).)'(^d 




^ '-'tot 


l + 2((^. 







(4.10) 



where contributions of 0{t^) to the total cross section are neglected. The factor -Fmt-iv 
shown in Table III is explained below. 

For the MT-III potential there is interference among the repulsive and attractive parts of 
the potential and the connection between the potential matrix element in forward direction 
and the total NN cross section is more involved. Since a^l^ = 0{-^) and the real part of 

the forward scattering amplitude t approaches a constant, both positive terms in Eq. (3.43) 
can again be combined as cr^/^ x F and we determined the factor multiplying as 
-^MT-III = 2.83, as shown in the last column of Table II. Because of the repulsion contained 
in the model MT-III the asymptotic limit is reached much later compared to the purely 
attractive potential model MT-IV. 

In both cases the rigorous asymptotic limits are 

^tot - ^.T X F, (4.11) 

where F is either -Fmt-iii or -Fmt-iv- The factors F contain the deuteron matrix element 
{fd\^\fd) and numerical as well as potential parameter constants. The factors F are larger 
than 2 for both potential models, and thus the rescattering process enhances the asymptotic 
cross section over the sum of the two NN cross sections, the shadowing has vanished already 
before that. 

We started from an expansion of the elastic n-|-d forward scattering amplitude ($|f7|$) 
in powers of the NN t-matrix. Due to the optical theorem the term first order in t ended up 
as a term second order in t in the total n+d cross section. The terms second order in t in the 
elastic forward scattering amplitude provide two terms in the total n+d cross section, one 
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of second order and one of fourth order in t. The one second order in t behaves hke O(^), 
whereas the one of fourth order in t decreases hke 0{\) in the hmit go going to infinity. It 
is therefore natural to group the terms together according to their power in t as we did in 
Eq. ( p.43| ), which coincides with their energy dependence. Thus the shadowing effect will 
disappear faster than the anti-shadowing effect. For the two model forces considered these 
anti-shadowing effects modify the naive expectation that the total cross section for n+d 
scattering tends to twice the total NN cross section. The true asymptotic result is larger 
than twice the NN cross section, which means that the two nucleons in the deuteron can 
never be considered to be independent. Finally it is obvious that the terms in Eq. (2.6), 
which are third order in t can not cancel the terms 0{t^) in the total n+d cross section. 
They may however modify the shadowing effect, since there will be also contributions of 
order O(t^) to the total n+d cross section. Thus the asymptotic value of the n+d total cross 
section in our model of spinless nucleons interacting by local forces is the first term on the 
right hand side of Eq. (4.9) in case of the purely attractive Yukawa potential MT-IV, and 
(j^^f X F with F = 2.83 for the potential model MT-III. These are exact results for our 
potential models. 

V. SUMMARY 

In view of new precise measurements of the total n+d cross section at energies above 
100 MeV nucleon laboratory energy and of precise solutions of the Faddeev equations using 
modern NN forces |jl| it is of interest to understand the leading rescattering effects which 
modify the naive expectation that the total n+d scattering cross section is the sum of the 
np and nn total cross sections. In a model of spinless nucleons we investigate the first 
few terms of the multiple scattering series for the forward elastic n+d scattering amplitude 
resulting from the Faddeev equations in the high energy limit. Although the treatment 
is purely nonrelativistic and enters far into the region, where relativity is important, we 
think it is interesting to know the asymptotic behavior at high energies for pure potential 
models without absorption. Absorption processes (particle productions) occurring in a rel- 
ativistic context will however change the results presented here 0,^. In accordance to the 
naive expectation we extract from the second order term in the multiple scattering series 
a shadowing effect proportional to the total NN cross section squared, which is negative 
and reduces the total n+d cross section. However, we also find a positive term proportional 
to the square of the real part of the NN forward scattering amplitude, which decreases in 
energy only proportional to the total NN cross section, whereas the shadowing term vanishes 
faster. Both terms are proportional to the expectation value of l/(nucleon distance)^ with 
respect to the deuteron wave function. 

In a numerical illustration using NN force models of Yukawa type the positive term has 
as limit the total NN cross section multiplied by a number larger than 2 (which would be 
the naive expectation). For the purely attractive potential this factor is 8.15, whereas for 
the potential with the additional repulsion this factor is 2.83. These are exact results for 
our choice of potential models. For both potentials the negative shadowing effect decreases 
faster as function of energy than the positive terms. 

We expect that the asymptotic expressions given in Eq. (3.43) start to be valid around 
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300-600 MeV depending on the potential model employed. Strictly speaking this can only 
be assured, if the Faddeev equation, Eq. (2.2), is solved exactly. The estimate given in this 
work is based on the exact evaluation of the term first order in the NN t-matrix in the 
n+d elastic forward scattering amplitude. From this calculation we found the corresponding 
asymptotic limit starting to be valid at the quoted energies. 

The formal expression, Eq. (3.43), has also been found in the context of the Glauber 
approximation 0,^]. A simple geometrical explanation of the deuteron matrix element has 
been given in Ref. 0. Our derivation is different and follows simply from the elastic forward 
n+d scattering amplitude evaluated in leading orders in the NN t-matrix and taking the 
high energy limit analytically. A similar investigation including spin and isospin degrees of 
freedom is in preparation. 
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APPENDIX A: DERIVATION OF THE MOMENTUM SPACE FORM OF THE 

INTEGRAL I 

In this Appendix we derive the limiting momentum space form of the integral I given 
in Eq. (3.20). Using the leading expression for the denominator as given in Eq. (3.29) the 
expression in Eq. (3.20) can be rewritten as 

3go J J Qz + q'^ + ie Att 

oo oo 11 /\//\ 

-1-1 

The X and x' integrals are elementary and one ends up directly with the expression given in 
Eq. (3.33). 
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TABLES 



TABLE I. Parameters of the Malfliet-Tjon type potentials. As conversion factor we use units 



such that nc= 197.3286 MeVfm=l. 





Va 


^iA [MeV] 


Vr 


m [MeV] 


{Vd\}iWd) [fm"^] 


MT-IV 


0.3303 


124.91 






0.6343 


MTTII 


3.1769 


305.86 


7.291 


613.69 


0.3090 



TABLE II. Contributions of the different leading order terms of the MT-III potential to the 
total n+d cross section as function of the NN laboratory energy. The last column shows the factor 
Fmt-III) which multiplies a^^^ and approaches a constant for higher energies. 



Eiab [MeV] 


2 ai^f [fm2] 


0(t2) [fm2] 


0(t4) [fm2] 


-''MT-III ~ „NN 


50 


63.945 


22.845 


-25.218 


2.71 


100 


23.438 


19.410 


-3.388 


3.66 


300 


9.885 


8.376 


-0.603 


3.69 


500 


7.314 


5.027 


-0.330 


3.37 


1000 


4.577 


2.595 


-0.129 


3.13 


1500 


3.400 


1.769 


-0.071 


3.04 


2000 


2.725 


1.341 


-0.046 


2.98 


2500 


2.283 


1.078 


-0.032 


2.94 


3000 


1.970 


0.899 


-0.024 


2.91 


3500 


1.738 


0.771 


-0.019 


2.88 


4000 


1.557 


0.673 


-0.015 


2.86 


4500 


1.414 


0.597 


-0.012 


2.84 


5000 


1.296 


0.536 


-0.010 


2.83 
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TABLE III. Contributions of the different leading order terms of the MT-IV potential to the 
total n+d cross section as function of the NN laboratory energy. The explicit expression for R is 
given in Eq. ( 4.1C| ). The definition for Fmt-iv is the same as in Table ||. 



Eiab [MeV] 


2 ai^f [fm2] 


0(^2) [fm2] 


0{t^) [fm2] 


-Pmt-iv 


R 


50 


82.057 


101.812 


-84.392 


4.48 


0.54 


100 


39.035 


76.478 


-19.098 


5.92 


0.71 


300 


12.076 


31.616 


-1.828 


7.24 


0.87 


500 


7.040 


19.590 


-0.621 


7.57 


0.91 


1000 


3.411 


9.989 


-0.146 


7.86 


0.94 


1500 


2.242 


6.693 


-0.063 


7.97 


0.96 


2000 


1.669 


5.035 


-0.035 


8.03 


0.96 


2500 


1.329 


4.035 


-0.022 


8.07 


0.97 


3000 


1.097 


3.349 


-0.015 


8.08 


0.97 


3500 


0.946 


2.894 


-0.011 


8.12 


0.97 


4000 


0.817 


2.513 


-0.008 


8.15 


0.98 


4500 


0.729 


2.240 


-0.007 


8.15 


0.98 


5000 


0.661 


2.034 


-0.006 


8.15 


0.98 
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FIGURES 



FIG. 1. The absolute values of the bound state wave functions ipd{k) as calculated from the 
two potential models MT-III and MT-IV given in Table 1. 

FIG. 2. The contribution of the first order term of the multiple scattering series for the n+d 
total cross section calculated exactly from the MT-III potential model as function of the asymptotic 
momentum |go (solid line). The dashed line represents the corresponding contribution of the first 
order term in the high energy limit. The variable |(?o = \J^^iab is chosen such that the kinetic 
laboratory energies for the two and three body systems are the same. 

FIG. 3. Same as Fig. 2, except that the two-body potential employed is the MT-IV potential. 

FIG. 4. Contributions to the n-|-d total cross sections as given by the MT-III potential as 
function of the corresponding NN laboratory energy. The solid line shows twice the total NN 
cross section. Successively added to this is the positive contribution (antishadowing) of 0{t^) (long 
dashed) and the negative contribution (shadowing) of 0{t'^) (dash-dotted). The magnitudes of 
those contributions are shown separately as dotted line for the positive and as short dashed line 
for the negative term. 

FIG. 5. Same as Fig. 4 except that the two-body potential employed is the MT-IV potential. 
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